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I discuss double-diffractive (double-elastic) production of the 77' and rj c mesons in the 
pp — > pXp reaction within the formalism of unintegrated gluon distribution functions 
(UGDF). The contribution of 7*7* — * rf fusion is estimated. The distributions in the 
Feynman xp (or rapidity), transferred four-momenta squared between initial and final 
protons [ti, <2) and azimuthal angle difference between outgoing protons (3>) are cal- 
culated and discussed. The results are compared with the WA102 data. Predictions at 
higher energies are presented. 



1. Introduction 

Recently the exclusive production of rf meson in proton-proton collisions was in- 
tensively studied close to its production threshold at the COSY ring at KFA Julich 
^and at Saclay^l. Here the dominant production mechanism is exchange of several 
mesons (so-called meson exchange currents) and reaction via Six resonance^. 

I present results of a recent study ™ (done in collaboration with R. Pasechnik 
and O. Teryaev) of the same exclusive channel but at much larger energies (W > 10 
GeV). Here diffractive mechanism may be expected to be the dominant process. In 
Ref!^ the Regge- inspired pomeron-pomeron fusion was considered as the dominant 
mechanism of the 77' production. Here I present results obtained in the formalism 
with unintegrated gluon distribution functions. Similar formalism was used recently 
to calculate cross section for exclusive Higgs boson productionl^El. There is a chance 
that the formalism used for Higgs can be tested quantitatively for exclusive (heavy) 
meson production where the corresponding cross section is expected to be much 
bigger. 

In Fig^ I show a sketch of the QCD mechanism of diffractive double-elastic 
production of rf meson (left diagram). For completeness, we include also photon- 
photon fusion mechanism (right diagram). 
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Fig. 1. The sketch of the bare QCD and photon-photon fusion mechanisms. The kine- 
matical variables are shown in addition. 

2. Formalism 

Following the formalism for the diffractive double-elastic Prod uction of the Higgs 
boson developed by Kaidalov, Khoze, Martin and Ryskin I^EH (KKMR) we write 
the bare QCD amplitude for the process sketched in Fig.l as 

-^pplp^p =*t / d k^ t V(k 1 ,k 2 ,P M ) , 2 .2 p 

J ""0,* ft l,t ft 2,t 

The bare amplitude above is subjected to absorption corrections which depend on 
collision energy. The vertex function V(k\, k 2 , Pu) in the expression Q describes 
the coupling of two virtual gluons to the pseudoscalar meson. The details concerning 
the function V{k\, k 2 , Pm) can be found in El 

The objects fgf/(xx, x[, fcg )t , kf it , h) and f°f/(x 2 , x' 2 , fc^, k 2it , t 2 ) appearing in 
formula Q are skewed (or off-diagonal) unintegrated gluon distributions. They are 
non-diagonal both in x and k\ space. Usual off-diagonal gluon distributions are 
non-diagonal only in x. In the limit xi j2 — ► x'i 2 , k$ t — > fc^ 2 t and £i, 2 — > they 
become usual UGDFs. In the general case we do not know off-diagonal UGDFs very 
well. It seems reasonable, at least in the first approximation, to take 

/^/(zr^UoVM,*!) = ^ f^i^kh) ■ f^ixukl) ■ F.ih) , (2) 
f° g f/(x 2 ,x' 2 , kl t ,kl,tM = y/f^ (4^o%) -/f (X2,kl t ) ■ Ftih) , (3) 

where Fi(ti) and F\{t 2 ) are usual Dirac isoscalar nucleon form factors and t\ and t 2 
are total four-momentum transfers in the first and second proton line, respectively. 
The above prescription is a bit arbitrary. It provides, however, an interpolation 
between different x and k t values. 

Neglecting spin-flipping contributions the average matrix element squared for 
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the p(7*)p(7*) — > Pfw/ process can be written as El 

I^C^p I 2 * ^ |F rr ^^(fc;, fc 2 2 )| 2 |k M | 2 lk 2 , t |W(cE) ■ (4) 

3. Results 

We have shown in RefF" that it is very difficult to describe the only exsisting high- 
energy (W ~ 30 GeV) data measured by the WA102 collaboration ^ in terms of 
the unintegrated gluon distributions. First of all, rather large cross section has been 
measured experimentally. Using prescription JJJ and on-diagonal UGDFs from the 
literature we get much smaller cross sections. Secondly, the calculated dependence 
on the azimuthal angle between the outgoing protons is highly distorted from the 
sin $ distribution, whereas the measured one is almost a perfect sin This 
signals that a rather different mechanism plays the dominant role at this energy. 




1 10 10 2 10 3 10 4 
W (GeV) 



Fig. 2. otot as a function of center of mass energy for different UGDFs. The 7*7* fusion 
contribution is shown by the dash-dotted (red) line. The world experimental data are 
shown for reference. 

In FigO I show energy dependence of the total (integrated over kinematical 
variables) cross section for the exclusive reaction pp — > prj'p for different UGDFs. 
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Quite different results are obtained for different UGDFs. This demonstrates huge 
sensitivity to the choice of UGDF. The cross section with the Kharzeev-Levin type 
distribution (based on the idea of gluon saturation) gives the cross section which is 
small and almost idependent of beam energy. In contrast, the BFKL distribution 
leads to strong energy dependence. The sensitivity to the transverse momenta of 
initial gluons can be seen by comparison of the two solid lines calculated with the 
Gaussian UGDF with different smearing parameter cto = 0.2 and 0.5 GeV. The 
contribution of the 7*7* fusion mechanism (red dash-dotted line) is fairly small 
and only slowly energy dependent. While the QED contribution can be reliably 
calculated, the QCD contribution cannot be at present fully controlled. 



Fig. 3. Two-dimensional distribution in t\ x ti for the diffractive QCD mechanism (left 
panel), calculated with the KL UGDF, and the 7*7* fusion (right panel) at the Tevatron 
energy W = 1960 GeV. 



In FigOU I present two-dimensional maps t\ x ti of the cross section for the 
QCD mechanism (KL UGDF) and the QED mechanism (Dirac terms only) for 
the Tevatron energy W = 1960 GeV. If |*i|, \t 2 \ > 0.5 GeV 2 the QED mechanism 
is clearly negligible. However, at |i2| < 0.2 GeV 2 the QED mechanism may 
become equally important or even dominant. In addition, it may interfere with the 
QCD mechanism. 

In Table 1 I have collected cross sections (in nb) for rf and r/ c mesons for W = 
1960 GeV, integrated over broad range of kinematical variables specified in the table 
caption. The cross sections for r\ c are very similar to corresponding cross sections 
for rj production and in some cases even bigger. 
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UGDF 


rf 


Vc 


KL 
GBW 
BFKL 
Gauss (0.2) 
Gauss (0.5) 


0.4858(+0) 
0.1034(+3) 
0.2188(+4) 
0.2964(+6) 
0.3793(+3) 


0.7392(+0) 
0.2039(+3) 
0.1618(+4) 
0.3519(+8) 
0.4417(+6) 


7*7* 


0.3095(+0) 


0.4493(+0) 



4. Conclusions 

The existing models of UGDFs predict cross section much smaller than the one 
obtained by the WA102 collaboration at the center-of-mass energy W = 29.1 GeV. 
This may signal presence of subleading reggeons at this "low" energy. 

Due to a nonlocality of the loop integral our model leads to sizeable deviations 
from the sin 2 $ dependence (predicted in the models of one-step fusion of two vector 
objects). 

The diffractive QCD mechanism and the photon-photon fusion lead to quite 
different pattern in the (tijia) space. 

Finally we have presented results for exclusive double elastic r\ c production. 
Similar cross sections as for rj' production were obtained. Also in this case the 
results depend strongly on the choice of UGDF. 

Measurements of the reaction(s) in the title would help to limit or even pin 
down the UGDFs in the nonperturbative region of small gluon transverse momenta 
where these objects cannot be obtained as a solution of any perturbative evolution 
equation, but must be rather modelled. 
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